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1. Introduction 

1.1 The survival probability 

The goal of this paper is to calculate the survival probability, taking into account the contribution of hard 
processes described by the BFKL pomeron exchange. The diffractive Higgs production is a typical hard 
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process, in which the Higgs is produced from the one parton shower due to gluon fusion. This process can 
be calculated in perturbative QCD. 



The signature of this process is the existence of so called large rapidity gaps (LRG), in which no par- 
ticles are produced (see Refs.|, |]). For the LHC energies and for diffractive Higgs production at the cm 
rapidity equal to zero, there are two rapidity gaps. The first is between the right moving final protons and 
the Higgs boson, the second is between the left fast moving proton and the Higgs boson. 



As was noticed by Bjorken M, in hadron hadron collisions, there is a considerable probability that 
more than one parton shower can be produced. Therefore, one needs to suppress such a multi parton 
shower production, since it can produce particles that fill up the rapidity gap. This suppression can be 
characterized by the survival probability || [| . 



To illustrate what survival probability is, it is instructive to calculate it in the simple eikonal model 
for soft pomerons. Soft pomeron means that there are no perturbative contributions from short distances, 
and only soft non - perturbative processes contribute to high energy asymptotic behavior. The survival 
probability is defined in the eikonal formalism as [Q, ^] 

, { o2^ _ f\M(8,b)\*e-W)<eb 
1 1 J\M(s,b)\W 1 j 

where M is the amplitude for the hard process under consideration, in impact parameter space (where b is 
the impact parameter), at the centre of mass energy yfs. In this paper, this is the amplitude of diffractive 
Higgs production from one parton shower. e~ n ^ gives the probability that additional inelastic scattering 
will not occur between the two partons at impact parameter b. £1 (b) is called the opacity or optical density. 
Therefore, the numerator is the amplitude for the exclusive process, while the denominator is the same 
process, due to the exchange of one pomeron. 



The survival probability was estimated in Ref.Q, in the eikonal approach for exclusive central diffrac- 
tive production at the LHC. The survival probability here was given for the process illustrated in Fig. |], 
in terms of the impact parameter b. Generally, in all these models the survival probability is given by the 
expression 



, | , _ I ^i rf2 6 2 {A H (h)A H (b2) ( 1 - Mih + b 2 )) 2 )) 2 ( | 



Aff(b) is the hard pomeron amplitude in impact parameter space b shown in Fig. |l[ The amplitude Ah (b) 
for hard pomeron exchange can be calculated in perturbative QCD, and is responsible for the production 
of two gluon jets, with BFKL ladder gluons between them (see Fig. [l]). In this model the hard pomeron 
in Fig. |l] emits the Higgs. An{b) is given in the impact parameter b representation by the expression |^] 

- ir~* (L3> 
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Figure 1: Central diffractive production in the two channel cikonal model in proton scattering due to pomeron 
exchange. 

R? H = 7.2GeV~ 2 . As shown in Fig. |], Au(bi) and Aufo) denote the hard pomeron amplitude above 
and below the Higgs signal respectively. The contribution A s , shown in Fig. [I], denotes the soft pomeron 
amplitude. ( 1 — A s ) includes all possible initial state interactions due to the exchange and interaction of 
soft Pomerons. ( 1 — A s ) also includes the possibility that the two initial nucleons in Fig. [j], do not interact 
at all. 

The survival probability was found to be 5 - 6 % for the single channel model. In the two channel 
model, the survival probability here is 2.7% at the LHC energy of yfs = 14000 GeV . The upper bound for 
the survival probability, in the constituent quark model (CQM), was found to be 6.0% ± 0.1% at the LHC 
energy. This is almost the same as the survival probability found in the single channel model. The two 
upper bounds, intercept at an energy just above the typical LHC energy. This suggests, that the upper 
bound for the survival probability, should be 2% — 3% for measurements at the LHC. 

The first attempt to estimate the contribution of hard (semi - hard) processes to the value of the 
survival probability, was made by Bartels, Bondarenko, Kutak and Motyka in Q. They considered the 
contribution of this "fan" pomeron diagram, to the value of the survival probability, and found that this 
contribution is rather large. Namely, the value ranges from 3.17% for a s = 0.15, to 1.6% for a s = 0.25, 
(where a s is the QCD coupling). 

The aim of this paper, is to calculate the BFKL pomeron (see Fig. ||), and the first enhanced diagram 
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for the BFKL pomeron (see Fig. |6|). These calculations are in the symmetric QCD dipole approach. Be- 
cause in proton proton scattering, there is no reason to assume that the mean field approximation, based 
on the "fan" diagram, can work. The ratio of the two contributions of Fig. || and Fig. || are calculated, 
and used to estimate the value of the survival probability. 



This paper is organised in the following way. In section ^, the coupling of the BFKL pomeron to the 
colour dipole (section |2.1| ), and the triple pomeron vertex (section ^^ ), are calculated in the momentum 
representation. Using these results, the BFKL pomeron amplitude shown in Fig. || is calculated (sec- 
tion ||]3|) , and the first enhanced diagram shown in Fig. ||, is calculated (section . 



Section || is devoted to the survival probability, estimated in the QCD dipole approach. The ratio of 
the two contributions of Fig. || and Fig. || is calculated, which is used to estimate the value of the survival 
probability (section |3.2| ). It turns out, that this ratio is not small, and indicates the importance of taking 



into account all enhanced diagrams. Therefore, in section all enhanced diagrams are summed in the 
toy model (see Ref.Q). The fact that the two dipoles have different sizes is neglected. From the calculation 
of the ratio of Fig. ||] and Fig. [6], the value of the parameter d of this model is determined, (d is the low 
energy amplitude for one pomeron exchange). Using this parameter, the value of the survival probability 
was estimated as the ratio of the diffractive Higgs production in this model, and Higgs production in one 
parton shower (for single pomeron exchange). It turns out that the survival probability is rather small. 



In the conclusion, the results for the value of the survival probability is presented. A discussion is 
given on the dependence of the value of the survival probability, on the choice of intercept of the BFKL 
pomeron. The significance of higher order hard rescattering contributions to the survival probability, is 
also discussed. 



2. The conformal eigenfunctions of the vertex operator and the triple pomeron vertex 

All calculations are carried out in the momentum representation, and the strategy and notation of Ref. || 
is closely followed. Firstly, the pomeron coupling to the QCD colour dipole is introduced (see Fig. |2|) in the 
momentum representation. Secondly, an explicit expression for the triple pomeron vertex (see Fig. ||) is 
derived. Using both these formulae, the BFKL pomeron (Fig. Q), and the first enhanced diagram (Fig. ||), 
are calculated in the symmetric QCD dipole approach. 

2.1 The BFKL pomeron vertex function 

The vertex coupling the BFKL pomeron to the couple dipole is illustrated in Fig. ^. Here, q is the 
momentum transferred along the pomeron, and x\2 is the transverse size of the dipole. In the notation of 
Refs. 0, [H , the eigenfunctions for the vertex in coordinate space are defined as 
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1 . . a 1 



(^01,^20 = -1 — (2.1.1 

\ X\2 J V ^12 / 



x n \ E^- U {xi 2 ) where Xij = Xi — Xj and Xi are the transverse coordinates. The conformal 

dimensions are defined as 

Figure 2: The vertex 
of interaction of Pomeron 
with the dipole. 

n n 
1= --iv 1 = --- %v (2.1.2) 

n is the conformal spin, and is an integer. The energy levels of the pomeron are the BFKL eigenvalues 
given by 

u(n,u) = a sX {l) =a s (2^(1) -^( 7 ) -^(l- 7 )) (2.1-3) 
where in this paper the notation 

a s = £ (2.1.4) 

7T 

is used, and where V'(Z) = d\nT (f) / df and T(/) is the Euler gamma function. Since the only intercept 
uj(n = 0, u) is positive at high energies, the contribution with can be neglected. Lipatov in Ref. |7| 

introduces the following mixed representation of the vertex. 



E^ u {x) = / d 2 Re^ R E n ' v (R+ |, |) (2.1.5) 



where 



_ 2^7H>r(3-«/ + a)r(§+«/) 
^-i-^r(f + ^ + i)r(|-^) (2 - L6) 



A more convenient expression of Eq. ( 2.1. 5| ) for the vertex was calculated in Ref. (8| and is given as 



£„"'" M = (OT*)'"2- 6 "T 2 (1-^) X (2.1.7) 
j t (^) * (^) " (-D " (^) ^ (qXh 

where J 7 are the Bessel functions of the first kind. In Eq. ( |2.1.7| ), q and q* are the components of the 
momentum q transferred along the pomeron, in the complex representation. That is 
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\Q\ = QQ 

where q = q x + iq y q* = q x - iq y (2.1.8) 

In order to work in the momentum representation when calculating the single pomeron amplitude (Fig. |j), 
and the first enhanced diagram (Fig. |6|), it is necessary to express the vertex function explicitly in the 
momentum representation. In Ref.|| it was shown that in the momentum representation, the vertex 
function is given by the following Fourier transform 

E « ,1 7) = / 4 «P I ^ exp BT M (2.1.9) 



X 



Here p denotes the momentum which is the conjugate variable of the dipole size x±2- The complex repre- 
sentation, to express the vector pm terms of its complex components p and p* (see Eq. J2.1.8 )), is used in 
Eq. (gXp . In Ref. @, this integral is written in the following factorised form 



E(p,qn) = h ^{q 2 )- iv 2-^Y\l-iu) ( E(p, q; j)E(p*, g*; 7 ) - E(p, q; -j)E(p*, q*; - 7 ))(2.1.10) 
where 



E(p*,q*;j) 



f dx ( q*x 

J 7^ {— 



E(p,q;j) 



dx* 



qx 



(2.1.11) 



At this point, it is assumed that n = 0, and hence 7 = 7 = — iv (see Eq. ( 2.1.2 )). This is because 
the only intercept u> (n = 0, v) is positive at high energies (see Eq. (|2.1.3[ )), so the contribution n / is 
neglected from now onwards. Let E(p* ,(/*; 7 ) and E(p,q;~/) be denoted as E(p*,q*;v) and E(p,q;v) for 
n = 0. After integration over x and x* , the expressions for E (q,p; v) and E (q*,p*;v) are found to be ||] 



E(p,q;u) 



E(j>*,q*;v) 
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(2.1.12) 
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Hence, using Eq. ( 2.1.12Q and the expression for 6 n =o,^ m Eq. ( p. 1.6 ), the RHS of Eq. ( 2.1.10 ) can be 
written in the explicit form as 



E (p, q; v) 



2 2 + 2i »(p 2 ) iv -h^T 2 {\-iv) Y{iu) 

v r 2 [ \ + w) r ( -iv) 



(2.1.13) 



v*"\ 2 



,113 1 q z \ „ / 1 1 3 1 

X 2-?] W i il/ , 1 — Zl/ , =■ 2-^1 ^ i W -, \ — IV , ; r-^ 

1 4 2 '4 2 ' '4p 2 y V 4 2 4 2 4(p*) 2 



+ 



r 2 (i + ^) r(-ii/) 



, , 1 1 ■ 3 1 . , . ^ \ _ / 1 1.3 1 . 1 . iq , 

X 2-^1 1 W, 1 IV , 1 + IV , - O-Tl 1 IV, 1 IV , 1 + IV , ; — - 

1 4 2 '42 ' ' 4p 2 y V 4 2 '42 ' ' 4 (p*) 2 



..*\ 2 



where E" (p, <f; n = 0, 7) is written as E (p, q; v). For single pomeron exchange, (see Fig. 0), the conformal 



spin v has opposite signs at the two vertices at the ends of the pomeron. In section 2.2 when Fig. H is 



calculated, it is assumed that q = to simplify the calculation. Hence, from Eq. (2.1.13), the product of 
two vertices, (for q = 0) takes the form 



1 1671"^ 

E(p,q = 0; v ) E(p,q = 0; -v) = — — — 

z/ z p L 



(2.1.14) 



2.2 The triple pomeron vertex 



q,k v 




q-q\,k,Vi Qi,k,Vi 

Figure 3: The triple Pomeron vertex. 




In this subsection the triple pomeron vertex, illustrated 
in Fig. [3| is calculated explicitly in the momentum represen- 
tation. It is defined in Refs. [U^, |ll| as an integral over 
the centre of mass position vectors ( xoi, xq2, X03) , and the 
conformal dimensions (7,71,71) as 



G3P [q,k,n = 0,7,71,72 ) = G 3P [q,k,v,vi,v 2 



(2.2.1) 



d 2 xi d 2 x 20 d 2 x 30 

Eg' {X 10 ,X20) E k ' {X20,X 30 ) E q l k {X30,X 10 ) 



^12^23^31 



To calculate the triple pomeron vertex explicitly, the mixed representation of Lipatov in Ref.|?J is 
used for the vertex eigenfunctions Eq ' u (see Eq. fl2.1.5p ). Note that to simplify the calculation of the first 
enhanced diagram of Fig. y, it is assumed that q = for the momentum transferred along the pomeron, 
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above and below the pomeron loop. Hence, the triple pomeron vertex shown in Fig. |3| is calculated for 
q = 0. In Ref.[[ll|] this mixed representation was used in the definition of Eq. (2.2.1) to give the expression 



G 3 p(q = 0,k,u,vi,u 2 ) 



1 

2^ 



^01 -2w-l„ ^ 



01 



01 



d x 2 



b 01 



yi 2 ,.yi 2 

x 12 x 02 



C 02 ^fc l X 02 J»12 



^_fc (,^12J 



(2.2.2) 

It is also assumed that v\ = zv 2 = in evaluating the expression of Eq. (12.2.2| ) for the triple pomeron 
vertex, for the following reason. When evaluating the integrals over v\ and u 2 , for the expression of Fig. [| 
(see Eq. (|2X2|) ), one expands the BFKL functions uj (v\) and u ( v?) around the saddle point v\ = v 2 = 
(see Eq. (|A-2-9|) ), which gives the largest contribution to the integration. In the appendix, the integral 
of Eq. ( p. 2. 2 ) is evaluated to give the triple pomeron vertex as an explicit expression in the momentum 
representation in Eq. (A- 1-14) as 



1 2~ 2iu T 3 ( - - iu) T 2 ( iu) 
G 3P (q = 0,k,u,u l ^ 0,^^ 0,) = fc*"-!- V 2 l . V ; (2.2.3) 



2.3 The single pomeron amplitude 




Figure 4: Central diffractive production in colour dipole scattering due to single pomeron exchange. 

In this subsection the single pomeron amplitude with Higgs production shown in Fig. ||| is calculated. 
The two dipoles are separated by a rapidity gap Y, and they have transverse sizes t\ = x\ 2 and r 2 = x' 12 . 
The momenta conjugate to the dipole sizes are pi and p 2 , and q is the momentum transferred along the 
pomeron. For simplicity it is assumed that q = 0. The single pomeron amplitude with Higgs production, 
in the QCD dipole approach is denoted Mniggs ( n = 1 , Y) , where Y is the rapidity gap between the two 
incoming protons, and n = 1 denotes the single pomeron exchanged between the two protons. In this 
notation, the single pomeron amplitude, with Higgs production of Fig. f|, has the expression 0, ||, 12, 13] 
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M Higgs (n = l,Y) = P^ KL (p 1 ,p 2 ,Y,q = 0)A H (5Y H ) 
where P BFKL (pi,p 2 , Y,q = 0) is the single pomeron amplitude given by the expression 

P BFKL (puP2,Y,q = 0) = ^ J ±-V(u)e^E (p 1 ,q = 0,u)E(p 2 ,q = 0,-u) 



(2.3.1) 



(2.3.2) 



and where Au(dYjj) denotes the subprocess contribution which produces the Higgs, such as the 
quark triangle subprocess of Fig. || The typical rapidity window, which the Higgs Boson occupies is 
5 Yjj = In I ^— ^ I , where m is the mass of the proton. The simplest subprocess with the largest contribution 



for Higgs production in the standard model, is the quark triangle shown in Fig. After the subprocess 
amplitude of Fig. |5| is contracted with the gluon propagators, the expression for the contribution of the 
quark triangle shown in Fig. || is given by 0, 15] 




Higgs 



Y u 



quark propagator 



Figure 5: Quark triangle subprocess for Higgs production. 



A H (SY H ) = A ( Mjj) [kfk 2 



where the factor A ( M\ ) has the value |T|, 0, |T|, |19| 



(2.3.3) 



A {Ml) 



2 l-a s {M 2 H ) (V2G F )' 



TT 



(2.3.4) 
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where Gf is the Fermi coupling. D (y) appearing in the single pomeron amplitude of Eq. ( |2.3.2| ) is 
given by 



V M = ( 2 V ln 2 ( 2 - 3 - 5 ) 

I" +4) 



(n = 0, v) is the solution to the BFKL equation defined in Eq. ( [2.1.3 ), where in the high energy limit 



one takes n = 0. From now on the notation u(n = 0, v) = u ( v) is used. Assuming that the conjugate 
momenta p\ and P2 of the two scattering dipoles in Fig. || are equal in magnitude, Eq. ( |2.1.14| ) can be used 



for the product of the two pomeron vertices. Hence, Eq. (2.3.5) can be written as 



P BFKL (P1 = P2 = P, Y,q = 0) = ^f ^—^— 2 e«W (2.3.6) 



p 2 J 2vri(i +l ,2) 

The integration over v can be evaluated at the saddle point v = of u (u). In this way, the RHS of 
Eq. ( |2.3.6| ) becomes 

p»«(ft = a = * r, ?= o) = ( j^m) * ^=" ,r 

Hence, the final expression of Eq. ( p.3.1| ) for the process of Fig. ^ reads 

M Higgs (n = l,Y)= ^ ( ^ { ^ =m ) ' e"^ ¥ A H (5Y H ) (2.3.8) 
This is the expression for the single pomeron amplitude, including Higgs production, of Fig. ||. However, 



Eq. ( 2.3.8] ) is written in the approximation that s = (p± + P2) 3> Mfj ■ Since we expect the Higgs mass 



to be large, we take into account the main correction due to this mass, namely, we make the following 
replacement. 

e ^=o)Y s (j_yW* _ |_Lj BHY ^p y s ^(y-HMl/^) (23Q) 

where m is the mass of proton., x\ and %2 are equal to k 2 /s\ and k 2 /s2 (see Fig. ||) with s\ = (pi + k\) 2 
and S2 = (^2 + &2) 2 )- Using the well known kinematic relation s\ S2 = Mfj s and since k 2 = k 2 = Mfj/2 



(see Ref. |J7j for example). Finally, Eq. (2.3.8) looks as follows 

^(^U)=^( ( ,„ (t/ ! 0)y) )' e ^=°Hr-HM^ m * ))AH{SYH) (2310) 
As one can see in Eq. (2.3.10) the single Pomeron exchange does not depend on the value of Higgs boson 



rapidity (Yh) but depends on SYjj = ln(M^/4m 2 ) which characterizes the window in rapidity occupied by 
the heavy Higgs boson. 
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2.4 The first enhanced amplitude 



In this subsection the amplitude for the first enhanced amplitude, with Higgs production shown in Fig. [6|, is 
calculated. The pomeron loop is between the two rapidity values Y\ and Y 2 . Hence, one needs to integrate 
over these two rapidity values. There is also an integral to evaluate, over the unknown momentum k in 
the pomeron loop. The enhanced diagram with Higgs production, in the QCD dipole approach is denoted 
Mniggs ( ti = 2 , Y) , where n = 2 denotes the splitting of the exchanged pomeron, into two branches forming 
the loop in Fig. ||. The amplitude of Fig. [f| is the first hard rescattering correction, to the single pomeron 
amplitude of Fig. ||. In this notation, first enhanced amplitude, with Higgs production of Fig. || is given by 



M mggs (n = 2,Y) = 2PZl L ced (Pi,P2,Y,q) A H (5Y H ) (2.4.1) 

where P^hancediPi'PZiY, q) is the BFKL pomeron amplitude for the first enhanced one - loop diagram, 
which has the expression given below in Eq. fl2.4.2| ). The factor of 2 in Eq. ( |2.4.1|) , 

comes from adding the 

two identical contributions of Fig. ^, due to the two ways the Higgs is emitted from the two branches of the 
pomeron loop. In order to obtain the complete contribution of Fig. ^, both possibilities for Higgs production 
from the two branches of the loop must be considered separately, and added. P^hanced(Pi'P2> ^ Q) * s gi ven 
by the expression (see Refs. || f2|) 



PenhancediPl Y, q) = B / dv dv ' dv x dv 2 d 2 k I dY l 

!y h +Uy h Jo 



Y 



Yh-^SYh 



dYo 



x E(p l7 q,u) V{v) e^ Y - Y ^G 3P (v,v 1 ,v 2 ,q = 0,k) 

xV(ux) V{u 2 ) e^^)+ w ^))( yi - y2 )-^ 2 ) ln ( M «/ 4m2 )G 3 p(-^, l / 1 ,^,g = 0,fc) 
x e ^')Y2 E (p 2 ,q,-is') (2.4.2) 



where B 



aiir 4 ( a, N. 



8 V 2tt 2 



(2.4.3) 



where Yh is the rapidity of the Higgs boson. Yh here is considered equal to zero in the cm. frame, 
restricting ourselves to the production of the Higgs boson at rest in the cm. frame since it is the most 
likely experimental kinematic, and 5Yh = ln(Mj^/4m 2 ) characterizes the rapidity window, occupied by 
the Higgs boson. 



Using the same assumptions as section 2.2, the integral over k in Eq. (2.4.2) is evaluated in the appendix 



with the result given in Eq. (A-l-17) as 



"1,^2- 



d 2 kG 3P (u,u 1 -► 0,i/ 2 -> 0,q = 0, k)G 3P (-u', v x -► 0,v 2 -> 0,q = 0,k) 

>Q 2 2(iu'-iu) 6 ^ _ ^ T 3 ( 1 + iv <^ F 3 ( 1 _ iv } p2 ( _ iv ,j T 2 ( {v) 



8 \v\\v\ 



w 



(2.4.4) 
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D( y>) e «>(y')Y2 j 






2 ; 9, vj 


r„ r~ 


— W- 



Figure 6: Central diffractive production in colour dipole scattering due to pomeron exchange with a hard rescattering 
correction. 



Note that the v\ and v 2 in the denominator of Eq. ( |2,4.4| ) cancel with T> ( v\) and T> ( i^) in Eq. (2.4.2) 
(see Eq. fl2.3.5| )). When inserting Eq. ( [2.4.4 ) into the right hand side of Eq. ( |2.4.2| ), the delta function 
allows the integration over v' to be evaluated to give the result 



P, 



BFKL 

enhanced 



(p 1 ,p 2 ,Y,q = 0) 



Bit 2 



dv dv\dv2 



Y H +\5Y B 



dY 1 







2 



dY 2 

)r-(i + 



IV 



IV 



ivfAv^v 2 sin 2 (ivir) 
x V{v 1 )V{v 2 )e {ui ^ +ul{v2 ^ Yl ~ Y ^ - ^ 1 < m hI^ 2 ) E (p 2 ,q,-v) (2.4.5) 

Now the integrals over v,v±,v 2 and the two rapidity values Y% and Y 2 need to be evaluated. Y\ and Y 2 
are the upper and lower rapidity values for the pomeron loop in Fig. |6[ The details of the integrations are 
given in the appendix in Eq. (A-2-3|)-Eq. ( A-2-11), and the final expression is given in Eq. (A-2-11) as 



PenhancediPl >P2iY,q— 0) — 

32j?vr 8 5Y H f (2u(v=0)) 4 (2u"(v=0)f \ M *=o) (V-f m( < ) ) 
(2a s ) 5 p 2 ^(v =0) ^ Y y 2 J e 
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where the constant B is given in Eq. ( 2.4.3 ). Therefore, the full expression for the diagram of Fig. || 
given by Eq. ( [2.4.1 ) takes the form 



M Higgs (n = 2,Y) = (2.4.7) 
64£vr 8 5Y H ( (2c [v = 0)) 4 (%f(v = 0)) 3 ^ M*=o) (y-\ U ) 

3. The survival probability in diffractive Higgs production in colour dipole scattering 
due to pomeron exchange 

3.1 The definition of survival probability 

In this section the survival probability of large rapidity gaps, in diffractive Higgs production is calculated, 
using the ratio of Fig. ^| and Fig. |6] (see and Eq. ( 2.3. 8| ) and Eq. ( |2.4.7| )). To guarantee that there will 



still be a large rapidity gap (LRG) between the protons after scattering, all hard rescattering corrections, 
that could give terms filling up the LRG, must be taken into account. The survival probability, is the 
probability to just have the exclusive Higgs production shown in Fig. and not to have any higher order 
hard rescattering corrections, such as the first enhanced diagram of Fig. ^. In other words , the survival 
probability is the ratio of the calculated cross section for the Higgs boson production to the one form the 
one pomeron exchange. Hence, the survival probability of the LRG, is calculated by subtracting the sum 
over all hard rescattering amplitudes from the single pomeron amplitude of Fig. ||, and dividing the result 
by the single pomeron amplitude of Fig. [| itself, to obtain the correctly normalised survival probability. 
Therefore, the survival probability is defined as 

, ig2| >= M Higgs (n=l,Y) -Zn=2(-l) n Mm 99 s(n,Y) 

' ' M Higgs (n = l,Y) { " ' ] 

where Mjjiggs ( n , Y) is the n th order hard rescattering correction. For example, in the case of n = 2, 
the first hard rescattering correction Mniggs ( n = 2 , Y) is the contribution of the first enhanced diagram 
of Fig. |6|, which has 2 pomeron branches, forming the pomeron loop. In general, M^iggs (n,Y) is the 
contribution given by the diagram which has n pomeron branches. In calculating the survival probability, 
if only the first enhanced diagram is taken into account, and corrections of the order n = 3 and higher are 
ignored, then the formula of Eq. ( |3.1.1| ) reduces to 



<\S 2 \ > 



M Higgs {n = 1 , Y) - M Higgs {n = 2 , Y) 
M Higgs (n = l,Y) 



M H 



n 



2, Y) 



M H 



n 



(3.1.2) 



The ratio 



,(n=2,Y) 



(n—i ,Y) ^ S calculated in the next subsection, in the symmetric QCD dipole approach 
(see Eq. (|3.2,4j ) in section 3.2 ). It turns out that this ratio is not small and, therefore, all enhanced diagrams 
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need to be taken into account. Using the toy model suggested by Mueller in Ref.Q, all enhanced diagrams 
are taken into account in the Mueller - Patel- Salam - Iancu (MPSI) approach (see Refs.[pO|, 21, The 
formula for the scattering amplitude in this model was suggested by Kovchegov in Ref. |23]| . 



3.2 The QCD dipole approach 

The survival probability of large rapidity gaps, in diffractive Higgs production in the QCD dipole approach, 
is the probability for the exclusive Higgs production of Fig. f|, with a large rapidity gap between the 
Higgs signal and the two emerging dipoles. To calculate the survival probability, all hard rescattering 
corrections which could fill up the large rapidity gaps must be subtracted from the single BFKL Higgs 
amplitude Mjjiggs ( n = 1 , Y) , and the result must be divided by Mm ggs (n = 1 , Y) . If only the first 
enhanced rescattering correction M Higgs ( n = 2 , Y) is taken into account, then the survival probability in 
the symmetric QCD dipole approach is estimated as 



|5 2, = M Higgs (n = 1 , Y) - M Higgs (n = 2, Y) 
<l l> M Higgs (n = l,Y) 



_ M Higgs ( n = 2 , Y) 

Mmggs ( n = 1 , Y) 



(3.2.1) 



where the amplitudes M(n = 1, Y) and M (n = 2 , Y) have been calculated in Eq. ( [2.3.8 ) and 



Eq. ( [2.4.7 ), respectively. Using the results of Eq. fl2.3.8|) and Eq. ( 2.4.71 ), then the ration j^|" = ^ Sp \ 



appearing in Eq. ( 3.2.1|) is found to have the expression 



M Higgs (n = 2, Y) 



Mmggs ( n 
2 Bit 5 



hY) 
5Y H 



(3.2.2) 



«2 (2a s fu"{v 







( 2u ( v = o) y 

Y 



2uj(v 



Y 2 




u»(v =o)Y y (U=0)Y 

2tt 1 ' 



where the constant B is given in Eq. ( 2.4.3j ). Here a typical value for a s , which depends on the mass 
of the Z particle, is used. It is expected that the Higgs will be produced with a mass of approximately 
100 GeV, which would give a value for the strong coupling constant a s ~ 0.12. This corresponds to a Z 
particle mass ||, of M z = 90.8 ± 0.6 Gev. 

The following values are to be found in Ref. [24| and Ref. [p5f |, for the strong coupling and the BFKL 
function. 



a s = 0.12 



uj{y = 0) = a s 41n2 



-y(is = 0) = Ua s ((3) 



C(3) = 1.202 



(3.2.3) 



Assuming that the rapidity gap Y at the LHC is 19, and using the numerical values given in Eq. ( |3.2.3 ), 
the right hand side of Eq. ( [3.2.2 ) then yields following. 
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M(n = l,Y) v ; 

This value is not small and increases with energy. Therefore, it shows that all enhanced diagrams have 
to be taken into account. In the next section all enhanced diagrams are summed in the toy model. 

3.3 The toy model approach 

In this subsection, the survival probability is calculated taking into account all enhanced diagrams. The 
toy model proposed by Mueller in Ref. ||, is a model for describing pomeron exchange in onium - onium 
scattering. In the toy model, the dipole wave function of an onium is described by the generating functional 
for dipoles (§) Z(Y, [«]): 

Z(Y - Y ; [u]) = (3.3.1) 

/n 
P n (Y; ri,&i, r 2 ,b 2 , ■ ■ ■ ,ri,bi, . . . ,r n ,b n ) u{r u bi) d 2 r^d 2 h 

71 = 1 1 = 1 

Here, P n are the probabilities to find dipoles with sizes and impact parameters bi at rapidity Y and 
u(xoi,b) is an arbitrary function of the dipole of transverse size xoi, at impact parameter b . In the toy 
model which we are going to consider here we neglect the dependence of u on the size of dipoles and their 
impact parameters (see Ref. ||) . In this model Z(xoi,b,Y,u) degenerates to the generating function and 
obeys the following evolution equation (see Ref. [0]) 

dZ{ d y U) =AZ 2 (Y,u)-AZ (Y, u) (3.3.2) 



where A is the pomeron intercept. In section |2.3| and section |2.4| , the pomeron intercept can be taken 
to be the BFKL intercept A = u> ( v = 0) to provide a matching with the BFKL Pomeron calculus. The 
initial condition for Eq. ( |3.3.2| ) is given by 

Z(Y = 0,u) =u (3.3.3) 
The solution of the toy model Eq. ( 3.3.2 ), which satisfies the initial condition of Eq. ( |3.3.3D is ||, 23] 



Z(Y, u) = —— (3.3.4) 

u + (1 — u) r 

Eq. ( |3.3.4| ) gives the sum over all "fan" diagrams. To generalise this result to the sum over all essential 
enhanced diagrams, the MPSI approximation is used to sum over all diagrams, with pomeron loops larger 



than In Ref. [23], the forward scattering amplitude in the MPSI approximation was written and has the 



2 

form 

D(Y,d) = l-e X p(-d-^^z(j,u\z(j,v)\ l , ,,. , (3.3.o) 
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where d is the dipole amplitude (0 < d < 1) at low energy. Substituting for Z(Y,u), the right hand 
side of Eq. ( 3.3.4| ) in equation Eq. ( |3.3.5| ), yields the following expression for D(Y , d) (2^] 



D(Y, d) = (-l) n D{n, Y, d) = - ^ n\ ( -l) n d n e nAY (l-e A ^) " (3.3.6) 

n=l n=l 

At large rapidity values, one can make the approximation 1 — e A ~ ~ — e A ~ , such that Eq. ( 3.3.6[) 



can be re- written as 



D(Y, d) = -]T(-l) n D(n, Y, d) = -J2n\(-l) n d n e nAY (3.3.7) 



n=l n=l 



In Eq. (|3A7| ), the n th term is the amplitude for n - pomeron exchanges. Hence, equation Eq. (|3.3.7| ) is 



the sum over all hard rescattering correction amplitudes for pomeron exchange, in onium - onium scattering. 



This approach is used in Refs. |20|, |21|, |22|. To include Higgs production in the toy model, one has to replace 
one of the n dipole amplitudes, by the contribution Ah(SYh) from the subprocess for Higgs production. 
The leading subprocess, is the quark triangle shown in Fig. ||. Hence, for each of the terms, a factor of n is 
included, to account for the possibility that the Higgs can be produced from any of the n pomerons. After 



Higgs production is included in the toy model of Eq. ( 3.3.7 ), the resulting amplitude takes the form 



D mggs (Y ,d) = - ]T ( -If D Higgs {n,Y,d) 



n=l 
oo 



^-lfcT-W^ 1 ^ 'A H (SY H ) = -^-D(Y,d)A H (5Y H ) 



dd 

71=1 



(3.3.8) 



where D Higgs (n , Y , d) = d n ~ x nine 11 AY A H {5Y H ) = ^ D (n , Y , d) A H (SY„) 

(3.3.9) 

The notation Duiggs (n,Y,<£) refers to the toy model BFKL pomeron amplitude, including Higgs 
production, with n pomeron branches. This should not be confused with the notation Mffiggs (n, Y) , 
which refers to the equivalent n th order term in the symmetric QCD dipole approach. The n = 1 term in 
Eq. ( 3.3.Sj ), corresponds to the single pomeron amplitude of Fig. [I[ The n = 2 term in equation Eq. ( |3.3.8| ), 



corresponds to the first enhanced amplitude of Fig. |6], with the hard rescattering correction of the pomeron 
loop. In section section 3J., the survival probability was defined by the expression given in Eq. ( 13.1.11) . 



Hence, in the toy model approach, the survival probability takes the form 
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'.y.y* 



n = l,Y,d) - En=2 (-l) n Dm ggs (n,Y,d) 
DHiggs ( n = 1 , Y , d) 



(3.3.10) 



Inspection of Eq. ( 3.3.8 ) shows that the numerator on the RHS of Eq. (|3.3.10 ) can be rewritten as 



DHiggs (Y , d) = D H iggs ( n = 1 , Y , d) - )^ ( -1)™ D Higgs (n,Y , d) 



n=2 



(3.3.11) 



Hence, the toy model formula for the survival probability of Eq. ( 3.3. 10[ ) becomes 



DHiggs (Y , d) 
DHiggs ( n = 1 , Y , d) 



&D(n = l,Y,d) 

Air 



(3.3.12) 



Typically, the Higgs signal will occupy a rapidity window 5 Yjj = In j-^ Therefore, in the toy model, 
pomeron exchange between scattering dipoles separated by a rapidity gap of less than 5Yjj, should be 
excluded for Higgs production. Therefore, the toy model amplitude Mniggs (n,Y,d) should be divided 
by the scattering amplitude Mniggs ( n , 8 Yh , d) , which gives the scattering amplitude for dipoles separated 
by a rapidity gap less than dYjj. Taking this into account, Eq. (3.3.12) is modified to give the survival 
probability for diffractive Higgs production within the rapidity window SYh, as 



<|5 2 l > 



(<L D (Y,d)) /(<L D (5Y H ,d)) 



Y,d))/(f d D(n 



l,SY H ,d)) 



(3.3.13) 



In order to calculate the survival probability using the expression of Eq. ( 3.3.10j ) , the value of the 
parameter d, appearing in the expression for D (Y , d) , must be determined. To do so, it is useful to 



refer back to the calculation of section 13.21, where the ratio 



MHzggs(n 2 , Y) calculated, in the symmetric 

M Higg3 {n=l ,Y) ' J 



QCD dipole approach (see Eq. ( 3.2.4D ). In order for the toy model to be consistent with the QCD dipole 
approach, the ratio calculated in Eq. ( 3.2.4 ), should be the same in the toy model. Setting n = 1, Eq. fl3.3.9| ) 
gives for single pomeron amplitude in the toy model 



D mggs (n = l,Y,d)=e AY A H (SY H ) 



(3.3.14) 



Setting n = 2 in Eq. ( 3.3.S| ), the first enhanced amplitude in the toy model is given by the following 
expression 



D H iggs(n = 2, Y , d) = -Me 2AY A H ( 5Y H ) 



(3.3.15) 



Therefore, (using Eq. (|3.3.14| ) and Eq. ( 3.3. 15|) ), the following condition is imposed 
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M Higgs (n = 2, Y) 
M Higgs (n = l,Y) 



Duiggs ( n = 2 , Y , d) 
Dmggs ( n = 1 , Y , d) 



-Ad z e 



2 Jl AY 



A H (SY, 



H 



de AY Afj (SYh) 



-Me 



Ay 



(3.3.16) 



Substituting for 



M HiggB (n=2,Y) 



jY^-, — ( n=1 y) * ne resu lt of Eq. ( |3.2.4| ) on the LHS of Eq. ( 3.3.16 ), and setting the 
pomeron intercept equal to the BFKL intercept A = u ( v = 0) , to be consistent with the QCD dipole 
approach, enables one to calculate a value for d in the toy model. One finds 



d = 0.7 



(3.3.17) 



One can now proceed to calculate the survival probability, by taking into account all higher additional 
hard rescattering corrections, using the formula of Eq. ( 3.3. 12| ). From Eq. (3.3.7), the expression for 
D ( Y , d) can be written as 



OO OO POO 

D(Y , d) = -J2nl(-de AY ) n = " E / ^e"* ( - dte Ay )' 

n=l n=l ^° 



1 " / *T+ju*r < 3 ' 318 > 







After changing variables to u = d g a y + t , then the RHS reduces to 



D(Y d) =1 _ eXp (^) r duer 1= exp(^) / _J_X 
[ ' ' de AY J ky u de AY V de Ay y 



(3.3.19) 



If one notes that in general -4- T ( , x) = — — -, then substituting for D (Y , d) , the RHS of 
Eq. ( |3.3.19| ) in the formula of Eq. ( |3.3,13| ) , gives the following expression for the survival probability. 



/e— ^ ( exp(^)r(0,^)(l±^)-d 

V 7 \ eX P de A ^H J 1 ' de A ^« J de*»*H ) ~ 

( M 2 

The typical rapidity window 6Yh, which the Higgs signal is expected to occupy, is oYh = In ( 
where M\ ~ lOOGeF. The typical rapidity gap is expected to be Y = 19 for the LHC energy of yfs = 
lATeV . Setting the pomeron intercept equal to the BFKL intercept, A = u{y = 0) = a s 41n2 ~ 0.34 
(see Ref . |25| ) , the value for the survival probability from Eq. ( 3.3.20|) , is found to be 



< |5 2 | > = 0.004 (3.3.21) 
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This gives the survival probability as 0.4%. However, the larger survival probability is obtained by 
abandoning the BFKL intercept u {v = ) ~ 0.34, and replacing the intercept with that of the soft 
pomeron, A = a s = 0.12. In this case, the RHS of Eq. ( p. 3. 20] ) using the soft pomeron intercept, gives 
the following value for the survival probability, 



< \S Z \ > = 0.23 (3.3.22) 

Hence, using the value for the soft pomeron intercept A = a s = 0.12, the survival probability is found 
to be close to 22%. Alternatively, using a higher value for the strong coupling a s = 0.25, and replacing 
the intercept with that of an upper limit for the soft pomeron intercept, A = a s = 0.25, the survival 
probability from Eq. ( 3.3.20| ), is found to be 



< \S 2 \ > = 0.022 (3.3.23) 
Which is considerably lower and close to 2%. This value is close to the value estimated by the Tel 



Aviv group in Ref. Q, and the Durham group in Ref. 27]. The values found for the survival probability. 



which depends on the choice of intercept are summarised below. 

Pomeron intercept A Survival probability < | 5 2 | > 

BFKL intercept A = to ( v = ) ~ 0.34 0.004 

Soft pomeron intercept A = a s = 0.25 0.022 
Soft pomeron intercept A = a s = 0.12 0.23 

Therefore, from these results it is clear that the survival probability depends critically on the intercept 
chosen. More specifically, the survival probability, as a function of the intercept A is not monotonic. The 
survival probability increases, as the intercept A decreases in value. For large rapidity gaps Y, then from 
the formula of Eq. (3.3.20), the survival probability is approximately proportional to 



<\S 2 \ >oc — . I — .. (3.3.24) 

11 exp(2A ( Y - 5Y H )) 

The typical LHC value for the rapidity gap Y between scattering dipoles is Y = 19, and for the predicted 
Higgs mass of Mfj, the rapidity window occupied by the Higgs, is expected to be S Yjj = In ( |^ 



Hence, provided Y — 5Yh > 0, then the expression of Eq. Q3.3.24 ) explains why, the survival probability 
increases as the intercept A decreases. 

Based on these results, in the toy model, the hard rescattering contributions from higher n corrections, 
range from 0.4% up to around 22%. Hence, the corrections are substantial and need to be taken into ac- 
count when calculating the survival probability, d in the toy model takes the value in Eq. fl3.3.17| ) d = 0.7. 
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This is less than unity. By inspection of the summation in Eq. ( 3.3.8| ), one can see that d is large enough, 
such that the terms n = 3 and higher, will give significant corrections to the survival probability calculated 
in this paper. 



To summarise, it is found firstly that d is large, giving significant higher contributions. Secondly, these 
higher contributions need to be taken into account, when calculating the survival probability. 

4. Conclusion 

The main results of this paper are the following. 

1. The first calculation of the enhanced BFKL diagram for diffr active Higgs production. 

2. Estimates for the survival probability for the full set of enhanced diagrams using the simplified toy 
model. 

3. The results of this estimate for the survival probability, show that the value depends crucially on 
the coupling constant of QCD, and that the multi pomeron exchange gives a substantial contribution 
to the survival probability. 

It was found that in the most consistent result for the survival probability, the value is rather small, 
0.4%. In conclusion, this paper shows that hard processes give a substantial contribution in the calculation 
of the survival probability. This paper is the first step forward towards obtaining reliable estimates of the 
influences of hard processes at high energy. 
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A. Appendix 



A-l Calculation of the triple pomeron vertex 



In this section the triple pomeron vertex is calculated to give an explicit expression in the momentum 
representation. This will be useful for calculating the first enhanced diagram of Fig. || in section 2A. In 



the expression for Fig. ||, (see Eq. ( [2. 4. 2D ), the BFKL functions oj(vi) and uj (z^) are expanded around 
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the saddle points v\ = v 2 = 0. This gives the largest contribution to the integration (see Eq. ( |A-2-9 )). 



Hence, in this subsection the triple pomeron vertex is calculated, in the limiting case when v\ = v<± = 0. 
It is assumed at the start of the calculation that v\ and v 2 are small and finite, however at the end of the 
calculation v\ and vi are put equal to zero. The triple pomeron vertex shown in Fig. || was defined in 
section (see Eq. ( ^2l|) ). 

G3P [q,k,n = 0,7,71,72) = G 3P (q,k,v,vi,v^ (A-l-1) 

E' (xio,x 2 o) Ej (x 2 o,^3o) E 'it {Xzo,X X q) 

^12^23^31 q 



A useful expression, to be found in Ref. |TlJ] , was given in Eq. (2.2.2) in terms of the mixed represen- 
tation of the vertex function E? ,v (x) (see Eq. ( 2.1. 5| )) as 



G 3 p(g = 0,k,v,vx,v 2 ) = — / —3— x 01 e 2 Ui 2 2 2 x 02 E k \x m )x l2 (^12) 

(A-l-2) 



In Eq. (A-l-2), it is assumed that q in Fig. [3| is zero. This is because for the calculation of the first 



enhanced diagram in section |2.4| (see Fig. |6|), the momentum q transferred along the pomeron above and 
below the loop, is set to zero, to make the calculation simpler. In Fig. [], there are two triple pomeron 
vertices, at opposite ends of the pomeron loop. Here, the momentum k is the unknown momentum in the 
pomeron loop. Evaluating the integral over xoi in Eq. ( A-l-2| ) gives an expression where the dependence 
on the momentum k is explicit, namely 



G 3P (q = 0,^,^2) = 2 3 " 2 ^-^ (fc 2 ) ^ ( \ w ^ g») gM% } (A.L3) 

r ( 1 + zzv + ivi + 1V2) 

where g 3 p (7,71,72) is the multidimensional integral related to the triple BFKL pomeron interaction, 
j;iven by [[Tl]] 

«, P ( 7 7] 72 ) = [ — / — / tUL \ R -b! 127+271+272-4 (A _ 1A) 

9iP{l, 11,72) J | x+ |2-2 7l | x _|2-2 72 J \ R+ \»n\R_\»rx J \R' + \^2\R'_\2y 2 \ n - n -\ 
where in the notation of Ref. (ll| , 

1 . 1 . 1 

7 = 2 + w 



x + % 



£± 

2 



71 


= 2 + 


iv 1 


72 


= 2 +W2 






n 


R + 




X- 


= x - 


' 2 




R\ 


= R' 


X- 


R'_ 








+ T 




2 



(A-l-5) 
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Consider the part of the integration over R' in Eq. ( [A- 1-4 ), which takes the form 



I (B! 



d 2 R' 



{R' + f 12 {R'_f 12 



R--R-) 



/\ 27+271+272-4 



dR' 



i\ 7+71+72-2 



R' + y 2 (R'_y 2 

dR'* 
R';Y 2 (R'*Y 2 



R.-RJ 



, 7+71+72-2 

R*.-R\) (A-l-6) 



In Eq. ( A-l-6 ) the complex notation d 2 R' = dR' dR'* has been used. Evaluating the integrations over 
R' and R'* gives 



/ 



dR' 



R — R ') 7+71+72-2 



1 

7T2 



{R' + y 2 {R'_y 2 

r ( I - 71 - 72 - 71) (R+-(x+-x) y^+''-l 
r ( 2 - 7 - 71 - 72) 



1 

X_2 



1 1 



+ r J I 2' 2' 71 +72 + 7 9 



1 i? + — ( X + — x) 



X- 



+ Tr-sr ( 2 - 7 - 71 - 72 ) r ( 7 + 7l + 72 - ~ ) x.( 7+7i+72 ) - 2 



(A-l-7) 



Inspection of the right hand side of Eq. ( A-l-7| ) shows that one has a singularity at 7 = 71 + 72 + 1, 
in the limiting case when v\ = V2 = (see Eq. ( A-l-5|) ). In this case, T(2 — 7 — 71 — 72) tends to infinity, 
which means that the first term on the RHS of Eq. ( [A-l-7] ) vanishes and the second term gives the largest 
contribution. Therefore, in this limiting case 



dR' 



R' + y< 2 {r'-Y 2 



R- ~ R-) 



A 7+71 +72 -2 



- 2 r ( 2 - 7 - 71 - 72) r ( 7 + 71 + 72 - ^ I 



(7+71+72) -2 

Inserting the result of Eq. ( |A-l-8 ) back into the result of Eq. ( A-l-4| ) gives 



(A-l-8) 



93 P (7)71,72) 



dR, 



dxA 



x+-n 



,7+71+72-2 



Rf J + xl 1 {x + - n) 72 ( R + - x+) 72 
HK) {{x* + )-nY + ^ 2 - 2 



R*) 72 



d(x* + ) 



(x* + ) 7x((x + ) -n) 72 ((i?;) -(x* + )) 



72 



1 / 3 

x -r 2 ( 2 - 7 - 71 - 72) r 2 7 + 71 + 72 - - 

7T V I 



(A-l-9) 



Now, using the notation for 7, 71 and 72 defined in Eq. ( [A- 1-5 ), <73 P ( 7, 71, 72) becomes in the limit 
that v\ = V2 = 0, 
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93 P (7,71,72) = .lim 



x lim 

iV2— * 



dR, 



R 



dR* 



dx-i 



x + — n) 2 



R* 



VlV2 



dx* 



R+-X+Y- 
1 



-IV\ 



( x%) * +lU2 (x* + -n) ( R* + - x * + ) 



-IV2 



(A-l-10) 



It is instructive to leave v\ and v 2 as small but finite in the indices, and let them be driven to zero at 
the end of the calculation, to avoid divergent integrals. Now integrating over x + and x* + gives the following 
result 



93p (7,71,72) 



lT 2 (--iv)T 2 (iv) lim f^k-x 



(A-l-11) 



1 1 - + iv\, - + ivr, 1 + 4^i, — ) + R i & 2iUl TT 2 Fi ( - 

r(l + 4wi) V 2 2 V 2 



Mui, ~+ivi, 1, -R + I x 



lim / —7- — 

7 (i?7 s ) 1+2w2 

1 



x 



— — — ^ , — - — 2-ti — h ^^'2, — h iv 2 , 1 + 4zi^2, + R-l 7t 2-n oivo-, — Vivi-, 1, ri . 

r(i + 4^ 2 ) V 2 2 R* + J y+J \2 ' 2 ' ' +y 

In the limit that iv\,iv<i — > then the factor r ( I t^J^Ls^^l ) — > 7;- anc [ Eq. ( A- 1- 1 1| ) reduces to 



93p (7,71,72) 



lim 



dR 



+ 



ivi->0j R 



x lim 

if 2— * 



? l+2ii/i 



2'2' 'i? 4 



2*1 ( ^,1,^-] + ^i^ 1 2^1 ( ^,1,#H 



1 1 

2' 2 ; 



, 1+2^2 



x vrr 2 ( t -i u ) T 2 {iu) 



(A-l-12) 



Finally evaluating the integral over R + in Eq. ( |A-1-12|) gives the result for g% p (7,71,72) as 



03p (7,71,72) 



1 
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Substituting this result for g%p ( 7, 71, 72) of Eq. ( |A-1-13| ) into the expression of Eq. (A- 1-3), the triple 
pomeron vertex is given explicitly in the momentum representation, in the limit that v\ = v 2 = 0, by the 
expression, 
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To calculate the first enhanced amplitude of Fig. ^, there is an integration to be evaluated, of the two 
triple pomeron vertices at both ends of the loop, over the unknown momentum k (see Eq. ( 2.4.2) )), which 
takes the form 



J d 2 k G 3P (q = 0, k, v, v x ,v 2 )G 3P (q = 0, k, -v', v x , v 2 ) 
Inserting the result of Eq. flA-1-14 ) gives 
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Now to integrate over fe, it is useful to make the change of variable I = Ink. Then the right hand side 



of Eq. (A-l-16|) reduces to a delta function in v and v', to give the result 



d 2 k G 3P (q = 0, k, v, v x ,v 2 )G 3P (-v' , v x ,v 2 ) 
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A-2 Calculation of the first enhanced amplitude 
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Once the integral over the unknown momentum k in the pomeron loop in Fig. ^ has been evaluated, the 
first enhanced diagram with Higgs production can be calculated from Eq. ( p. 4. 2 ). Inserting the result of 
Eq. (|A-T17| ) in the right hand side of equation Eq. Q2.4.2 ), gives 
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where B= -^(^ " 
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The largest contribution to the integral over v in Eq. ( A-2-lj) , is when 
pole at iv — > |, such that 
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It is assumed, that the conjugate momenta p\ and p 2 of the two scattering dipoles in Fig. |6|, are equal. 
Using the expression of Eq. (2.1.14) for the two pomeron vertices, Eq. (A-2-1) reduces to 
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One has the singularity as iv — > \, in the integrand of Eq. ( A-2-3 ), due to the factor - — ^tte- To 
remove this singularity, the substitution 

( a s (Y-Y 1+ y 2 ) \ 
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is used. The integrand in Eq. ( |A-2-3j ) then simplifies to 
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There is still the singularity in the integrand due to the factor 1 2 . To remove this singularity, it 
is useful to change the variables such that 
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Then Eq. ( [A-2-5 ) simplifies to 
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Integrating over u, the right hand side is found to be proportional to a delta function in the rapidity, 
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The integration over the two rapidity variables Y\ and Y 2 is now simple, because of the delta function 
in the integrand. The BFKL functions u and u {v 2 ), can be expanded around the saddle points v\ = 
and v 2 = 0. From the definition of Eq. ( [2.3.5 ), T> {fi) and T>{v 2 ) would vanish at v\ = and v 2 = 0. 



However, the factors of v 2 and v\ are canceled by the v 2 and v 2 , appearing in the denominator in the 



integrand of Eq. ( A-2-8 ). Hence, the integrals over v\ and v 2 can be evaluated at the typical values v\ = 



and v 2 = 0, to give the result 
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Finally, there are the integrations over v\ and iaj to evaluate, which are just Gaussian integrals. 
Therefore, the result for the right hand side of Eq. ( |A-2-9| ) is 
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Taking the rapidity Y to be 19 for the LHC energy y/s = 14000 GeV , and u{v =0) = 4a s In 2, 
then the first and second terms in the brackets of Eq. ( |A-2-10| ) are the largest terms, and hence at leading 
order 
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